Limits on Universality in Ultracold Three-Boson Recombination 
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The recombination rate for three identical bosons has been calculated to test the limits of its 
universal behavior. It has been obtained for several different collision energies and scattering lengths 
a up to 10 5 a.u., giving rates that vary over 15 orders of magnitude. We find that universal behavior 
is limited to the threshold region characterized by E < h 2 /2fj,\2a 2 , where E is the total energy and 
fli2 is the two-body reduced mass. The analytically predicted infinite series of resonance peaks and 
interference minima is truncated to no more than three of each for typical experimental parameters. 

PACS numbers: 34.10.+x,32.80.Cy,05.30. Jp 
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The development of Bose-Einstein condensates (BECs) 
with tunable properties makes possible condensates with 
a wide range of interaction strengths. Several experi- 
ments P, UM El Hi have investigated these properties by 
changing the atomic interactions with an external mag- 
netic field near a diatomic Feshbach resonance. As is 
well-known, these interactions are characterized at low 
temperatures by the two-body scattering length a, which 
covers the full continuum of positive and negative values 
near the Feshbach resonance. In fact, all of the essen- 
tial properties of BECs are determined by the scattering 
length. 

Because the two-body scattering length is the only rele- 
vant parameter, the precise shape of the two-body poten- 
tial does not matter for the many-body physics. One of 
the remarkable results to emerge from recent work on ul- 
tracold collisions is that this property holds even for the 
low temperature three-body physics. This shape inde- 
pendence has allowed theorists to choose any convenient 
two-body potential that reproduces the desired scatter- 
ing length. We will again exploit this freedom in the 
present work to study three-body recombination. Three- 
body recombination is the process by which three free 
atoms collide to form a diatomic molecule and an un- 
bound atom, setting free enough kinetic energy to make 
both atom and molecule escape from typical traps. 

The universal behavior of the three-body recombina- 
tion rate makes it possible to derive analytical expres- 
sions for a > [1 Q, i| and for a < 0. In the former 
case, theory predicts an infinite number of minima in the 
rate as the scattering length goes to positive infinity; and 
in the latter case, an infinite number of maxima as the 
scattering length goes to negative infinity. The physics 
behind both features is closely related to the Efimov ef- 
fect I n fact, it has been suggested that measuring 
the recombination rate while tuning through a Feshbach 
resonance might make possible some of the first direct 
experimental evidence of this intriguing effect. 

Tuning through such Feshbach resonances can dramat- 
ically limit the density and lifetime of BECs, however, 
since the three-body recombination rate was predicted 
0000 E3 — and recently verified experimentally 



— to increase with the scattering length as a 4 . More re- 
cently, three-body recombination has been used to create 
composite bosons by pairing fermions in ultracold gases 
[eH Il3| . The ultimate goal of this endeavor has recently 
been achieved with the observation of Bose-condensed 
pairs of fermion atoms [l4j . Despite its importance and 
recent advances, much work remains for the theory of 
three-body recombination. In particular, the universal 
behavior of the recombination rate has not yet been 
tested by accurate calculations. 

In this Letter, we show that the recombination rate 
for identical bosons is universal only for collision ener- 
gies in the threshold regime. Generically, the thresh- 
old regime is characterized by fc|<z| < 1, or equivalently, 
when the collision energy is the smallest energy in the 
system. For positive scattering lengths, the energy scale 
is set by the two-body binding energy; for negative scat- 
tering lengths, by the height of a potential barrier or, 
in some cases, by a two-body shape resonance. There- 
fore, for a fixed total three-body energy E, the relation 
E < E\2 indicates when the system is in the thresh- 
old regime where universal behavior is expected. In this 
expression, E\i = Ti 2 /2^i2a 2 is the two-body binding en- 
ergy and /ii2 is the two-body reduced mass. 

The experimental consequences of restricting the range 
of universal behavior are striking. At any nonzero tem- 
perature, rather than observing an infinite series of reso- 
nances or minima, only a finite number of either will be 
observable as the scattering length is scanned from — oo 
to +oo — and even those will be washed out. For in- 
stance, at 1 nK, only three resonances and three minima 
can be hoped to be seen. 

We will further show that the analytical formulas de- 
rived in Refs. 0, S hold only at zero energy. At 
finite energies, when the scattering length is tuned out of 
the threshold regime (E > £12), the analytical formulas 
break down because they do not take proper account of 
three important finite energy effects: unitarity, thermal 
averaging, and higher partial waves. Unitarity limits the 
rate to finite values at finite temperatures for large scat- 
tering lengths, and leads to a saturation effect jH 0|. 
Thermal averaging takes account of the fact that exper- 
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iments are performed at fixed temperature rather than 
fixed collision energy, and higher partial waves must al- 
ways be included, in principle. A generalized Wigner 
threshold law guarantees that the J T = + contribu- 
tion dominates at threshold, where J is the total orbital 
angular momentum and 7r is the overall parity. The next 
leading contribution, 2 + , grows with energy as E 2 and 
with scattering length as a 8 , and can quickly become 
comparable to the + rate. 

We obtain the recombination rates by solving the 
Schrodinger equation numerically usingthe adiabatic hy- 
perspherical representation (see Refs. 0, lis lilt UM for 
details of our implementation) . The key to this approach 
is that the dynamics of the three-body system are re- 
duced to the motion on a set of coupled effective po- 
tentials that depend only on the hyperradius R. The 
hyperradius is a collective coordinate that represents, in 
some sense, the overall size of the system. The effective 
potentials are determined by solving the adiabatic equa- 
tion 

H^^n) = U u (R)^ u (R;Cl) 

where Q denotes the five hyperangles representing all de- 
grees of freedom besides R. The adiabatic Hamiltonian 
i? a d includes the kinetic energy for these hyperangles as 
well as all interactions. The effective potentials U V (R) 
are then used in the radial equations (atomic units will 
be used unless otherwise noted), 
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where F v is the hyperradial wave function, E is the total 
energy, and the three-body reduced mass \i is related to 
the atomic mass m by \i = m/y/3. The nonadiabatic cou- 
pling Wuv' is responsible for inelastic transitions such as 
three-body recombination. The effective potentials give 
a very intuitive picture for these complicated systems. 
Moreover, the calculations can be made as accurate as 
desired by including more channels in the equation above 
(all rates quoted here are accurate to at least three dig- 
its and were obtained with seven channels). The radial 
equations Q are solved using the variational i?-matrix 
method 19] in order to extract the S-matrix. 

The three-body recombination rate Ks is defined in 
terms of the 5-matrix as 0, E E3, E3| 
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where k = y/2/iE is the hyperradial wave number, and i 
and / label the initial and final channels, repectively. 
The present results were obtained using the mass of 
helium atoms and the model dimer potential v(r) = 
Dsech 2 (r /ro) with D and tq adjusted to give a single 
two-body s-wave bound state. 



In the adiabatic hyperspherical representation, recom- 
bination for a > is driven primarily by the broadly 
peaked nonadiabatic coupling between the lowest three- 
body entrance channel and the highest molecular chan- 
nel 0- In this picture, there are two indistinguish- 
able pathways for recombination, leading to the so-called 
"Stuckelberg oscillations" . This interference phenomena 
modifies the a 4 dependence of the rate, suppressing it for 
certain values of a. At zero energy, the analytic results 
predict these minima to be equally spaced on a loga- 
rithmic scale and separated by a factor of approximately 
e*! - w 22.7, where a = 1.0064. 

For a < 0, the recombination rate is enhanced for 
particular values of a and, with the help of the adia- 
batic hyperspherical representation, can be interpreted as 
three-body tunneling through a potential barrier in the 
entrance channel . The nonadiabatic coupling is local- 
ized at small R behind this barrier so that recombination 
is suppressed for energies below the barrier maximum. 
That is, unless the collision energy matches the energy 
of a three-body resonance trapped behind this barrier. 
Under these conditions, transmission through the barrier 
jumps and strong enhacement of the recombination rate 
can be observed. As with the interference minima, these 
resonances can be associated with Efimov physics 10]. 
and are also predict to be equally spaced on a logarith- 
mic scale (separated by a factor of about 22.7). Figure^] 
shows the J 71 " = + recombination rates calculated at en- 
ergies in the range 0.1 nK to 1 mK. Figure^a) shows the 
first three resonance peaks for a < 0, and Fig.QIb) shows 
the first three interference minima for a > 0. For small 
values of the scattering length, the rates for all energies 
lie along a common, universal curve. For any given en- 
ergy, the rates depart from this universal curve at some 
value of the scattering length, with the highest energies 
departing soonest. 

In the adiabatic hyperspherical representation, the an- 
alytic recombination expression is derived under the as- 
sumption that the collision energy is in the threshold 
regime. It is natural, then, to conclude that the ana- 
lytic expression — and thus the universal behavior — is 
only valid in the threshold regime. The collision energy 
is in the threshold regime when it is smaller than other 
characteristic energies of the system. One obvious energy 
scale is the two-body binding energy E 12 . 

In Fig. the vertical dashed lines mark the scat- 
tering lengths determined from the relation E < E±2 
for each energy. It is clear that the two-body bind- 
ing energy provides a reasonable estimate for the do- 
main of universal behavior, i.e., for each energy, the 
rate curve for a less than this limit follows the com- 
mon curve. For a < 0, a better, more restrictive limit 
can be determined from the adiabatic potential since the 
threshold regime in this case requires energies less than 
the potential barrier maximum, t/ max = 0.079//za 2 |ll| . 
which reduces the limiting a by about a factor of three. 
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FIG. 1: Three-body recombination rate for (a) a < and (b) a > 0. The horizontal lines represent the unitarity limit un for 
each energy (reported as a temperature using E = ksT). The vertical lines represent the boundary of the threshold regime ajv. 
The solid curve represents the analytical formula Eq. The insets show the third resonance peak (a) and third interference 
minima (b) from a — 0. 



For example, in a 23 Na condensate at a temperature of 
100 nK, the recombination rate is expected to be uni- 
versal only for -3200 a.u.< a < 8650 a.u.; for 87 Rb, for 
-1650 a.u.< a < 4450 a.u. All of these values are well 
within the range that are already experimentally accessed 
near Feshbach resonances. 

We also show in Fig.^the analytical results 0, 01 0> E| : 

( 4590 (a 4 /m)smh(2i),) „ 

— J sin 2 [a ln(3|a|/2r )+*+1.63]+smh 2 77, U V U ' (3) 

I 360 {a 1 Ira) sin 2 [a ln(3a/2r )+$] a > 

where $ and 77* are unknown parameters. $ repre- 
sents an unknown small- R phase H (related to A* 111 
Refs. l8j,|9() and is chosen to give the best fit to the third 
interference minimum at 0.1 nK. The additional 1.63 rad 
of a < phase is predicted in 0. The value 77* =0.1 was 
found to give the best fit for a < 0. There is generally 
very good agreement with the numerical results for large, 
positive a/ro, and Eq. © appears to be essentially exact 
for zero energy recombination. It relies on the effective 
range expansion, however, and gets increasingly worse as 
\a\ decreases due to order r /a errors (here, r = 15 a.u.). 
The agreement is more qualitative for a < due to the 
small shift of the resonance peak positions. We found, 
though, that a 15% change in the extra a < phase 
gives good agreement with the 0.1 nK curve. 

One factor left out of Eq. J3J is unitarity (although a 
"unitarized" version has been proposed [15j to help ex- 
plain the experimental results in Ref. 0]). As the collision 
energy grows large compared to E12 for a fixed scatter- 
ing length, the probability of recombination approaches 
unity for the + partial wave. More relevant for experi- 
ments, unit recombination probability is also reached as 
the scattering length is increased at fixed collision energy 
E. The horizontal dashed lines in Fig. ^ denote the uni- 



tarity limit — un = 1927r 2 //ifc^r, obtained from Eq. (0 
by setting |S'| 2 =1 — for each energy shown. From the 
figure, it is clear that the recombination rate reaches the 
unitarity limit for positive a outside the threshold regime. 
For negative a, however, while the rate does saturate, it 
does so at a value about a factor of ten below unitarity. 
The main effect of unitarity is to restrict the number of 
resonances or minima observable at a given energy. 

A second factor neglected in Eq. © is the thermal 
average. Experiments are performed at fixed tempera- 
tures rather than fixed energies, so the thermal average 
becomes crucial for proper comparison with experiment. 
In the threshold regime, the recombination rate is con- 
stant as a function of energy, so the thermal average has 
no effect. Since we consider exactly the situation when 
the system is no longer in the threshold regime, thermal 
averaging can have significant effects. The thermally av- 
eraged recombination rate is 

{Ki)[T) = 2{^Tf J K ^ E ) E ^' E ' kBTdE ^ ( 4 ) 

where ks is Boltzmann's constant. Figure [21 illustrates 
the effects of thermal averaging at 0.1 /iK and 1 /iK near 
the second resonance peak and second interference min- 
ima. For energies solidly within the threshold regime, 
thermal averaging has little effect. For energies on the 
border of the threshold regime, however, averaging re- 
duces the intensity of both the peaks and minima, mak- 
ing their observation much more difficult. 

A third factor not included in Eq. JSJ is the contri- 
bution from higher partial waves. The J n = 2 + rate 
was calculated for -3000 a.u.< a <8000 a.u. and ener- 
gies from 0.1 nK up to 1 mK. The 2 + threshold law is 
K% (x E 2 a 8 0|, so for a finite energy, there will be a scat- 
tering length for which the 2 + contribution is comparable 
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FIG. 2: Thermally averaged recombination rate (a) near the 
second resonance peak and (b) near the second interference 
minimum. Circles and diamonds represent K3 and (if 3), re- 
spectively, at 0.1 ^K (open symbols) and 1/iK (filled symbols). 

to + . Figure 13 shows the thermally averaged 2 + rate at 
10 nK and 1 fiK along with 0+ for a > 0. (For a < 0, 
the 2 + recombination rate is many orders of magnitude 
smaller than for + , making it completely negligible for 
the present range of scattering lengths and energies.) It is 
clear from the figure that the 2 + rate dominates + at the 
second interference minimum for 1 /iK so that the total 
rate will show just one minimum. At 10 nK, the 2 + rate 
is merely comparable to + at the second minimum, cut- 
ting its depth in the total rate. For energies below 10 nK, 
the 2 + recombination rate is negligible in this range of 
scattering length; for larger scattering lengths, however, 
the 2 + recombination rate can contribute substantially. 
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FIG. 3: Comparison of the and 2 recombination rate for 
a > and temperatures of (a) 10 nK and (b) 1 /iK. 

Taken together, the unitarity limit, thermal averag- 
ing, and higher partial waves restrict the analytic results 
for ultracold three-body recombination to the threshold 
regime, i. e., when E < E\2- 

Experimentally, the consequences are rather dramatic. 
If we imagine tuning the scattering length using a Fes- 
hbach resonance, then a will, for instance, change from 
its background value to +00, then to —00, then again 
to its background value — all while the system is at es- 
sentially the same temperature. The analytic expressions 
predict that the three-body recombination rate grows like 
a 4 as the resonance is approached, goes through an in- 
finite number of minima as a — > +00, then has an infi- 



nite number of resonances as a returns from —00. Each 
series of features reflects Efimov physics, so measuring 
them might reveal evidence for this effect. The present 
calculations show, however, that the infinite series are 
truncated to a small number (w a/ir ln(3a c /2ro), where 
a c = %j \Z2iiyiE\2) for typical experimental parameters 
and that the contrast of the surviving features may be 
considerably reduced. The recombination rate is thus 
not a good candidate for observing physics related to the 
Efimov effect except at extremely low temperatures. 

Even though we have shown that the universal behav- 
ior described by existing analytic expressions is limited 
to the threshold regime, scattering lengths up to a few 
thousand atomic units are included. Moreover, a new 
sort of universal behavior dictated by the unitarity limit 
may take over and modifications to the analytic expres- 
sions along these lines have already been proposed 
Since we have used only one model potential, we are not 
in a position to discuss any universal behavior outside of 
the threshold regime. We expect, however, that recombi- 
nation for a > will be much as we have shown in Fig. ^ 
since it takes place at large distances where differences in 
the two-body potential will have little effect. For a < 0, 
the situation is just the opposite since recombination is 
a small distance process. The resonance positions as well 
as the a —> —00 limiting rate will likely then depend on 
the two-body potential. 
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